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Abstract 

This paper is devoted to the classification of integrable Davey-Stewartson type equations. A 
list of potentially deformable dispersionless systems is obtained through the requirement that such 
systems must be generated by a polynomial dispersionless Lax pair. A perturbative approach based 
on the method of hydrodynamic reductions is employed to recover the integrable systems along with 
their Lax pairs. Some of the found systems seem to be new. 
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1 Introduction 



In this paper we continue the programme of classification of integrable dispersive 2 -|- 1-dimensional 
equations started in [H [2] and present the complete classification of integrable two-component systems 
of second order 

ut = F{u,v,w, Du, Dv, Dw) 

n( n n n ^' ^J^^""^- ^'^■^> 

Vt = Lr[u, V, w, L)u, Dv, Vw) 

Here u{x,y,t),v{x,y,t),w{x,y,t) are scalar variables, Du, Dv, Dw denote the collection of partial 
derivatives of u, v, w with respect to x, y up to the second order and F, G are polynomials in derivatives 
with coefficient depending only on u, v, w. 

The approach proposed in pm2j consists of two steps: 

• One first classifies integrable dispersionless equations which may potentially occur as dispersion- 
less limits of systems in consideration; 

• One then reconstructs the dispersive corrections preserving integrability. 



One of the most famous examples within class (jl.ll) was introduced by Zakharov in [13] 

Ut = {uv)x + Uxx, Vt = VVx +Wx - Vxx, Wy = Ux, (1.2) 

together with the associated Lax pair 

Ipxy + ^Vtl^y + ^Ulp = 0, + i^xx + ]^W1p = 0. (1.3) 

Two more known examples |10^ [TT] include 

Ut = WUx +Vx+ Uxx, Vt = {vw)x - Vxx, Wy = Ux 

and 

Ut = {uv)x + WUx+Uxx, Vt = {vw)x + VVx - Vxx, Wy = Ux 

and their Lax representations are given by 

ilp^y + V1p - luipx = 0, + i^xx - Wipx = 

and 

2lP^y + VlPy - UlPx =0, + IPXX - WIPx = 0. 

The classification results obtained in this paper include 6 more integrable systems of the form (11. ip 
(Theorem 3), some of which seem to be new. 

The dispersionless counterpart of system ()1.2p read as 

Ut = iuv)x, Vt = VVx + Wx, Wy = Ux (1.4) 

and the corresponding Lax representation is 

SxSy + ^vSy + = 0, St + Sl + ^w = 0. (1.5) 
These can be obtained from (jl.2p and (jl.3p through a change 

s 

X ^ ex, y ^ ey, t ^ et, ip = (1-6) 

and taking the limit e — ?> 0. 

We recall that the dispersionless system (jl.4p is integrable by the method of hydrodynamic reductions 
(see e.g. [3|), a method which can be applied to any first-order quasilinear system 

A{u)ut + B{u)ux + C{u)uy = Q, (1.7) 

where A,B,C are n x n square matrices and u is an n-dimensional vector. Integrability in this sense 
means that for any fixed number A^, an integrable system of the form (jl.7p possesses an infinite class 
of solutions of the form 

u = u{R^,...,R^), V = v{R\...,R^), w = w{R^,...,R^), (1.8) 

parametrized by N arbitrary functions of a single argument, where the Riemann invariants i?* (also 
called phases) are assumed to satisfy pairs of commuting diagonal systems 

4 = ^*(i?)iii, Ri = X\R)Rl i = l,...,N. (1.9) 

In order that relations (jl.Sp and (jl.Op represent a nontrivial solution, it is readily seen that the 
velocities A* , /x* must satisfy the so-called dispersion relation 

D{X\n^) = det{X'A + B + fi^C) =0, i = l,...,N. (1.10) 

Throughout this paper, we consider only systems for which (ll.lOp defines an irreducible curve. More- 
over, compatibility of systems (|1.9p requires that the velocities A*(i?), //*(i?) obey the conditions 

T-^ = — iy^j, (1.11) 
A* — A-? /i* — ^ 

where A*- = djijX^,iij = dj^jiJ^. We admit the following definition of integrability of dispersionless 
systems 



Definition 1. A quasilinear system (1.1) is said to be integrable if for any number of phases N it 
possesses infinitely many N -phase solutions parametrised by N arbitrary functions of one variable. 



Owing to the fact that the flows (jl.9p are automatically commuting for = 1, one-component 
reductions of system (|1.4p are parametrized by the following relations 

u = R, v = V{R), w = W{R), Rt = iRV' + V)R,, Ry = j^^R,, (1.12) 

where V{R) is an arbitrary function and W{R) is fixed up to an integration constant by the relation 
W' = R{V')'^. Moreover, the normalization u = R is assumed to hold through a Miura-type trans- 
formation. It was shown in [H [21 [3] that deformability of hydrodynamic reductions can be used as 
a definition of integrability of dispersive equations. Indeed, formulae (II. Sp . (jl.Op governing hydrody- 
namic reductions can be deformed through the addition of a dispersive correction of the form 

u = u{R\ i?"), V = V{R\ i?") + e(. . .) + e2(. ..) + ..., 

w = W{R\ i?") + e(. . .) + e\. ..) + ..., (1.13) 



Ri = X\R)Rl + 6(. . .) + e\. ..) + ..., Rl = /.^(i?)4 + e{. . .) + e\. . .) + 



where the terms at are homogeneous differential polynomials in the x-derivatives of R^ of total 
degree k + 1 with coefficients being functions of the Riemann invariants. Using this construction the 
following definition of integrability of 2 + 1-dimensional dispersive systems was introduced in [1] . 

Definition 2. A 2 + 1-dimensional system is said to be integrable if all hydrodynamic reductions of 
its dispersionless limi^ can be deformed into reductions of the corresponding dispersive counterpart. 

To accomplish the classification program, the polynomial systems (jl.ip are written in the quasilinear 
form 

ut = aux + fiuy + 71;^. + 5vy + pw^ + e(. . .) 

, , , ^ , ^ ^ I Wy=Ux, (1.14) 

vt = <pux + WUy + r]Vx + TVy + KWx + e(. . .) 



under which they are often referred to as Davey-Stewartson type systems. Here the greek symbols 
a, . . . , denote functions of u, v and the nonlocal variable w = Dy^D^u and the terms at e are assumed 
to be homogeneous differential polynomials of total degree 2. The dispersionless limit 

Ut = + l3Uy + ^Vx + 5Vy + 

-■■ - Ux, (1.15) 



Vt = (t)Ux + 1pUy + r]Vx + TVy + KWx 

is thus obtained by simply setting e = 0. To classify integrable systems ()1.14p we first have to classify 
integrable dispersionless systems (ll.lSp and then classify dispersive deformations of the obtained 
systems. Classification of integrable systems (jl.lSp can be performed using the conditions of existence 
of 2- and 3-component hydrodynamic reductions (Section 2). The list of such systems is very extensive 
and moreover most of examples are not deformable to ()1.14p . However we notice that if a dispersive 
system (I1.14p is integrable then it possesses a Lax representation 

LV^ = 0, {dt - A)i^ = 0, 

where L and A are some linear differential operators with coefficients being functions in u, v, w and 
their derivatives. The dispersionless limit (|1.6p brings the Lax representation to a form 



F{Sx,Sy) — 0, St — G{Sx,S, 



where F and G are polynomials in Sx, Sy with coefficients depending on u, f , w. Thus if a dispersionless 
system is deformable then it possesses a polynomial dispersionless Lax representation and it is enough 
to restrict ourself to only such dispersionless systems. We therefore adopt the following classification 
approach: 



^which is supposed to be linearly non-degenerate, see Remark[2l 



• We first classify integrable dispersionless systems whicli possess polynomial dispersionless Lax 
representations (Section 2). In this paper we restrict ourselves to second degree polynomials; 

• We then classify dispersive deformations of such systems (Section 3). 

We finally emphasize that the same perturbative approach can be used to obtain the Lax pair of an 
integrable dispersive system from that of its dispersionless limit, given that this limit is integrable in 
the sense of hydrodynamic reductions. In Section 3 we present a construction which allows to deform 
dispersionless Lax pairs into their dispersive counterparts and give the complete list of integrable 
systems (jl.l4p along with associated Lax representations. 

2 Classification of dispersionless systems. 

In this section we describe and apply the two approaches to the classification of integrable dispersionless 
systems 

( ut = aux + (3uy + -yv^ + 6vy + pw^: _ , , 

[ Vt = (pu^ + IpUy + r]Vx + TVy + KWx 

where we assume that coefficients a, . . . , n are functions u,v,w. 

On one hand, the method of hydrodynamic reductions provides a sufficiently restrictive integrability 
criterion which permits to distinguish integrable systems within the considered class. Similarly to the 
3-soliton criterion for solitonic systems, it can be seen that the presence of 3-phase solutions (solutions 
of the form (jl.Sp . (jl.9p with = 3) is necessary and sufficient to prove integrability of dispersionless 
systems. However, requiring directly the existence of 3-phase solutions for system (j2.ip is hardly 
workable since the dispersion relation of system (j2.ip defines a cubic curve which may not be a priori 
a rational curve. The following theorem allows to overcome this difficulty by providing necessary and 
sufficient conditions for the existence of 2-phase solutions for any system (jl.7p . 

Theorem 1 (Ferapontov, Khusnutdinova [7J). The Haantjes tensor of an arbitrary matrix Q = 
{kA + B)~^ (lA + C), /c, / G M, is zero if and only if the following hold. 

• System ^1. 7\ ) possesses double waves parametrized by four arbitrary functions of a single argu- 
ment. 

• The characteristic speeds A*, /i* of two-component reductions are not restricted by any algebraic 
relations other than the dispersion relation U.10\) . 

Recall that for a given square matrix Q = (ujj), the Haantjes tensor is defined as [8] 

H^, = Niy^u^l - - N^.u^Xj + ^>>P^ (2-2) 

where the Njj^ are the components of the Nijenhuis tensor 

TV], = u;]du.ul - uldu^u;) - ^(S,.^^ - d^.u;]). (2.3) 

Thus the vanishing of the Haantjes tensor imposes a large set of first order differential constraints 
on the ten unknown functions a, . . . ,k. These conditions are very restrictive and allow to obtain a 
wide list of systems (I1.15P which possess 2-component reductions. Imposing the further condition of 
existence of 3-component reductions one then obtains the complete list of integrable dispersionless 
systems (jl.lSp . The list is extremely lengthy and most of systems are not deformable to (I1.14p . 

Remark 1. More precisely, it was possible to solve completely the case p = 0, while the case p ^ 
seemed to be out of reach for the moment. Nevertheless, it is remarkable that all deformable integrable 
systems obtained through our second approach fall into the class p = 0. 



Remark 2. In order to apply the perturbative approach it is necessary to exclude linearly degenerate 
dispersionless systems as their dispersive counterparts do not inherit the hydrodynamic reductions. 
These systems are characterized by the fact that their A^-phase solutions are continuous and do not 
admit the gradient catastrophe. A different perturbative approach would be required. Fortunately, 
these systems can be easily discarded by verifying an algebraic criterion that was introduced in [5]. 



Now we proceed to the classification of dispersionless systems (|2.ip which possess polynomial disper- 
sionless Lax representations. It was noted in the Introduction that the system 

ut = {uv)x, vt = vv^ + Wx, Wy = Ux, (2.4) 

possesses a dispersionless Lax representation 

SxSy + ^vSy + \u = 0, St + Sl + ^w = 0. (2.5) 
Introducing notations p = Sx, q = Sy, s = St the relations (j2.5|) can be rewritten as 

F := pq + -vq + -u = 0, s = G := —p^ — -w (2-6) 
and the compatibility conditions Sxy = Syx , Sxt = Stx , Syt = Sty take the form 

Py = qx, Pt = Sx, qt = Sy. (2.7) 
It is easy to see that relations (|2.6|) and (|2.7|) can be written as 

Ft = {G,F} mod F, equations ([23]), (2.8) 



where 



SG F} — _ + dGd£ _ d£dG 

dp dx dp dx dq dy dq dy 



and ^ = Ux^ + Vx^ , ■ ■ ■, i.e p,q are treated as symbols in (12. Sp . 

We now introduce the following definition of polynomial dispersionless Lax pair (see eg \12>\ 114]). 

Definition 3. A pair F{p,q),G{p,q) of polynomials in p,q with coefficients depending on u,v and 
u,v,w correspondingly is called a generator of a dispersionless Lax pair for a system \2. 1\) if 



1. F,G satisfy 

Ft = {G, F} mod F, £ (2.9) 



2. rank -tt^ -jt^ = 2, where K, — „ ,■ „ ,- 1 

\ du /jj=i^...^degF ^ dp'dqJ\p=q=0 

A polynomial Lax pair for a given system is not unique. Apart from transformations p — )■ fj,ip+fi2, q — >■ 
fiiq + fi3, s — )• /iis + //4, Hi, . . . fi4 = const we also have that if F, G are generators of a dispersionless 
Lax pair for a system then f{F), G are also generators of a dispersionless Lax pair for the same system 
for an arbitrary polynomial /. Moreover G is defined modulo F. Condition 1 is the compatibility 
condition while condition 2 implies that the Lax pair is faithful, i.e. the compatibility condition is 
satisfied due to both equations in the system (j2.ip . 



We now present the list of dispersionless systems which possess polynomial dispersionless Lax pairs. 
We restrict the classification of polynomial Lax pairs to those of degree 2 and leave the consideration 
of systems with Lax pairs of higher degrees for future studies. 



Theorem 2. If the system ^2. 1\) possesses a polynomial dispersionless Lax pair with degF = 2, degG 
mod F = 2 and if F is an irreducible polynomial then up to the group of invertible transformations 
V — 7> /(n, v), X ^ X, y ^ y + sx, u ^ u,w ^ w — su, s = const, and rescaling of variables it is one of 
the following ones: 



ut 

vt 

Ut 
Vt 

Ut 
Vt 

Ut 
Vt 

Ut 
Vt 

Ut 
Vt 

Ut 
Vt 

Ut 
Vt 

Ut 
Vt 



■- [UV)x 
VVx + Wx ' 

WUx + Vx 
{vw)x ' 

: WUx + {uv)x 
{vw)x + VVx ' 

: WUx + Vx- ^VVy 



{vw)a 



hvv^ 



WUy + Vx + VVy 



■ UUx 
2vUx + v'^Uy + uv 



WUy 



' W^ 



[v — hu)ux + h{w — ^v'^)uy + uvx — huvv 



„ . „ y 
-2hvux + h'^v'^Uy + {v — hu)vx + h{w — ^v'^)vy + Wx ' 

: {hu + v)ux + h{w — 2huv — ^v'^)uy + uVx — hu{hu + v)vy 
—2hvUx + h'^v'^Uy + {v — hu)vx + h{w — ^v'^)vy + Wx ' 

: {uv)x + WUx — ^v'^{l + 2hu)uy — huv{l + hu)vy 
{vw)x + Wx — 2hv'^Ux + h'^v^Uy — ^hv'^{l — 2hu)vy — 2huvvx 

■ {w — 2huv)ux + h^uv'^Uy + (1 — hu'^)vx + hv{u^ — l)vy 
—2hv'^Ux + h'^v^Uy + {w — 2huv)vx + h^uv'^Vy + vwx ' 



(2.10) 
(2.11) 
(2.12) 
(2.13) 
(2.14) 
(2.15) 
(2.16) 
(2.17) 
(2.18) 



where Wy = Ux, h = const. 



By the deg G mod F we understand the minimal degree of G modulo F. Parameter h in the above 
list, if not zero, can be scaled to 1. Note that if /i ^ then systems (I2.15p . (I2.16P reduce to (12.10p . 
while systems (j2.17p . ()2.18p reduce to (j2.12p and (|2.1ip correspondingly. The corresponding generators 
for the dispersionless Lax pairs are given in Table [H 



Eq no 



pq + ^vq + lu 

pq - \up + \v 

pq - \up + \vq 

pq + \hvq + \v — \up — ^huv 

pq + ^uq + j+vq"^ 

hpq — ^(t; + hu)q — h^vq^ + 

pq + (7 — hvq^ + \u 

pq+^{l- hu) - hvq^ + ^ - \vq{l 

pq — hvq^ — ^up + | 



G 



(!2T0]) 
(12TB 
(I2T2]) 
(!2T3]) 
(I2T4I1 
(I2T5D 
(I2T6D 
(I2T7D 
(|2T8]) 



-p - ^w 
-p^ + wp 
-p^ + wp 

-P^ +Wp- Tah^v'^ 



-p^ + v'^q^ 



16' 
wq 



hu) 



-p^ + h'^v'^q^ + \h{2w + v'^)q 
-p^ + h'^v'^q^ - {h^uv -hw + \hv'^)p - \w 
-p^ + h'^q^v'^ — ^huv"^ — ^hq{2hu — l)w^ + wp 
-p^ + h'^v'^q^ + wp- \hv'^ 



Table 1: Generators for the Lax pairs of dispersionless systems (|2.10p - (j2.18p . 

We finally make the following observation which greatly simplifies the computations entering the 
perturbation scheme in the next section. 



Proposition 1. The dispersion relation I11.10\) of every system in 112.10\) - ^.18\) is a rational cubic 
curve. 



Indeed the dispersion relation (ll.lOp of systems of the form (I2.ip defines a homogeneous cubic polyno- 
mial in the variables A*, and can be viewed as a ternary cubic form by introducing triples X'^Y'^ Z* 
through the relations A* = X'^/Z^, /i* = /Z^. It can be shown that for each of systems (j2.10p - ()2.18p 
the dispersion relation is a singular cubic form and therefore admits a rational parametrization. This 
can be verified by noting that the discriminant of the curve vanishes identically. This discriminant 
takes the form of a sixth-order polynomial and there exists a straightforward algorithmic procedure 
originally developed by Hesse [9J to obtain its explicit expression (see also the appendix in [12]). The 
singular point can either be a double point or a cusp and it should be noted that for systems ()2.15p - 
(j2.18p . the nature of the singularity can change in the limit /i — )• 0. For these cases, taking the limit 
in the corresponding Lax pair generators does not lead to a Lax pair of the limiting system. 



3 Dispersive deformations. The classification theorem. 

In this section we study the dispersive deformations of systems (I2.10p - ()2.18p up to order e. We illustrate 
the deformation procedure by recovering Zakharov's system (|1.2p and its Lax pair (jl.Sp from their 
dispersionless limits. To this end, dispersive corrections in the form of infinite series are added to 
equations ()1.12p namely 

u = R, v = V{R) + evi{R)R^ + e\.. .) + ..., w = W{R) + ewi{R)R^ + e^{.. .) + ... , 

Rt = [RV + V)R, + e {aiiR)R,, + a2iR)Rl) + e'(-) + • • • > (3.1) 

Ry = J^^R- + ^ {Mm.. + Mml) + eH-) + ■■■ 

where W = R{V')'^ and the unknown functions vi{R),wi{R), ai{R), Ai{R), . . . are to be determined. 
Again, it is assumed that the relation u = R remains undeformed through a Miura-type transformation. 
Equations ()3.ip are then assumed to form a formal solution of some second-order dispersive system 

ut = {uv)x + e {aiUxx + . . . + cr22UyWx) , Vt = vv^ + + e i(J23Uxx + • • • + a^iUyWx) , (3.2) 

where the nonlocal relation Wy = Ux is kept undeformed and ui, . . . , CJ44 are functions of u, v, w to be 
determined. Noting that V{R) is an arbitrary function in (jS.ip . we now require that equations ()3.ip 
define a solution of (13. 2p for any function V{R). This requirement together with the compatibility 
condition Rty = Ryt specify uniquely all unknown functions, up to a rescaling of e, and Zakharov's 
system (jl.2p is obtained. Up to the first order, the formal series solution is given by 

u = R, v = V{R)-e{V'fRx + 0{e^), w = W{R) + 0{e^), 
R^ = {RV' + V)Rx + 0(e2), Ry = j^^R. + O(e'), 

while higher order terms can be determined algorithmically. We point out that a sufficient condition 
for defining completely the dispersive deformation, i.e. such that no function remains arbitrary, is that 
the dispersionless limit must not be linearly degenerate (see Remark[2]). Moreover, it turns out that all 
systems in the considered class possess a unique integrable deformation. We do not exclude however 
the possibility that a given dispersionless system may possess two or more inequivalent dispersive 
counterparts and our algorithm would lead to all such deformations. Applying the same procedure 
to systems (|2.10p - (j2.18p . we obtain the following exhaustive list of integrable Davey-Stewartson type 
systems. 



Theorem 3. Up to the group of point transformations u ^ u,v ^ f{u,v) and re-scaling of variables, 
the following systems constitute the complete list of dispersive deformations of systems I2.10\ )- f2.18\) 



up to order e. 



vt = vv^ + Wx- evxx ' 

Ut = wux + Vx + euxx 
Vt = {vw)x - eVxx ' 

Ut = {uv)x + WUx + eUxx 

Vt = {vw)x + Wx - eVxx ' 

Ut = wux - \h?vvy + Vx + e{uxx + hv^ 
Vt = {vw)x + hvvx - eVxx 



(3.3) 
(3.4) 
(3.5) 
(3.6) 



Ut = -wuy + uux + A{v) + e {A^{u) + UyA{v)) , . 

Vt = v'^Uy + 2vux — wvy + uvx + Wx — e {A'^{v) — VyA{v)) 



(3.8) 



Ut = {v — hu)ux + h{w — ^v'^)uy + uVx — huvvy + e [B^{u) — huyB{v)) , 
Vt = —2hvUx + h?v'^Uy + {v — hu)vx + h{w — |w^)'(^y +Wx — e [B'^{v) + hvyB{v)) 

Ut = {hu + v)ux + [hw — 2h?uv — \hv'^)uy + uVx — hu{hu + v)vy 

+e(i?^(n) — 2huB{vy) — 2huxVy — huyVx + Sh'^vUyVy + 2h?uVy), (3.9) 



Vt 



-2hvux + h?v'^Uy + {v — hu)vx + h{w — \v'^)vy + Wx — e {B'^{v) + hvyB{v)) 



Ut = {uv)x + wUx — \v'^{l + 2hu)uy — huv{l + hu)vy 

+e(i?^(n) — 2huB{vy) - 2huxVy — huyVx + 2>h?vUyVy + 2h'^uVy), (3.10) 
Vt = {vw)x + Wx — 2hv'^Ux + h?v'^Uy — ^/iu^(l — 2hu)vy — 2huvvx — e (-B^(u) + hvyB{v)) 

( Ut = {w — 2huv)ux + h?uv'^Uy + (1 — hu'^)vx + hv{hu'^ — l)vy + e i^B'^{u) — huyB{v)) , , 

\ Vt = -2hv'^Ux + h'^v^Uy + {w - 2huv)vx + h?uv'^Vy + vwx - e {B'^{v) + hvyB{v)) ^ ' ' 

where Wy = Ux and A = Dx + vDy, B = Dx — hvDy. Parameter h in the above list if non-zero may be 
rescaled to 1. 

Note that several of these systems have already appeared in other classification works. Zakharov's 
system (j3.3p was recovered in [Tlj along with systems (|3.4p and (j3.5p . Moreover the systems 



Ut = 2wUx+Uxx, Vt = 2wVx-Vxx, Wy = {u + v)x, (3-12) 

and 

Ut = 2wux + Uxxi Vt = 2wvx - Vxx, Wy = {uv)x, (3.13) 



also appeared in |10] and were shown to be equivalent to (13. 4p and (13. Sp through a Miura transformation 
in [TT]. To the best of our knowledge, systems ()3.6p - ()3.1ip never appeared in other classification 
problems and seem to be new. 

3.1 Deformation of dispersionless Lax pairs 

Lax pairs for systems ()3.3p - ()3.1ip are obtained simultaneously through an analogous procedure. 
The dispersionless Lax pairs for systems (|2.ip are given by a pair of equations 

F{p,q) = 0, s = G{p,q) (3.14) 

and by relations 

Py = Qx, Pt = Sx, qt = Sy. (3.15) 



We can view the dispersionless system (I2.ip together with (13.14113.15]) as a quasihnear system of 
hydrodynamic type of the form (jl.7p with a reducible dispersion relation (jl.lOp . This system possesses 
an infinite class of solutions of the form 



u 



u{R\...,R^), V = v{R\...,R^), w = w{R\...,R^) 

p = p{R\...,R^), q = q{R\...,R^), s = s{R\ . . . , R^) 
Ri = fi\R)Ri, Ri = X'{R)Ri, i = l,...,N. 
Similarly to the construction of deformations for the system (j2.ip and its hydrodynamic reductions 



we seek the deformation of the dispersionless Lax pair as 

p = p{R\...,R^)+eO+e'0 + --- , q = q{R\ . . . , R"") + eQ + e\) + ■ ■ ■ , 
s = siR\...,R'') + eO+e\) + --- , 

where at we have homogeneous differential polynomials in derivatives of R^ of degree k + 1. Defor- 
mations of equations (|3.14p are in the form 

F{p, q) + e{fipx + f2Py + fsQy + + fbUy + hv^ + fjVy) = 0, ,^ 
s = G{p, q) + e{gip^ + g2Py + gsqy + g4Ux + g^Uy + g^v^ + gjVy + gsw^c), 

where = fi{u,v,p,q), gi = gi{u,v,w,p,q) and relations (|3.15p are kept undeformed. Notice that we 
need to deform equations ()3.14p up to the order e only. 



Dispersive deformations of the system (j2.ip are already constructed and thus it is only necessary 
to find series (|3.16p satisfying (j3.15p and (|3.17p . Analogously to the procedure of constructing the 
dispersive corrections for the equation one can algorithmically construct series p.l6p and obtain a set 
of differential equations on the dispersive deformations fi,gi. 

Equations (|3.17p . (|3.15p can be viewed as a nonlinear form of a Lax representation for corresponding 
dispersive system. To obtain the linear Lax representation we recall that p = Sx, q = Sy, s = St,ip = e~ 
and thus we can apply the transformation 

p = eDxlog{i^), q = eDy logiiP), s = eAlog(V'), (3.18) 

to relations (I3.17|) . 

Finally, the algorithm of deformations of dispersionless Lax pairs is the following: 



• For a given dispersionless system one obtains its dispersive deformations; 

• One deforms the corresponding dispersionless Lax pairs as in ()3.16p . ()3.17p using the disper- 
sive deformation of the system and its hydrodynamic reductions to obtain the nonlinear Lax 
representation; 

• One applies the linearization transformation. 



In the following Tables [2] and [3] we present nonlinear and linear Lax representations for systems (|3.3p - 
(|3.1ip . The transformation between the tables is of the form (I3.18p . 

We note that for a given dispersive system the deformation of the corresponding dispersionless Lax pair 
may be not unique. Moreover for systems ()3.3p -( f3rTT]) there exists 3 different dispersive deformations 
of their Lax pairs. In the tables we present the simplest deformations only as we need only one Lax 
representation for every system (I3.3p - (l3.1ip . 



Eq no 


Nonlinear pair 


(13:31) 


+ ^vq + {u + eq^ = 0, 
s+p^ + \w + epx=Q 


dm) 


pq - \up +\v + eqx = 0, 
s + — + epx = 


(|33]) 


- ^up + ifg' + eqx = 0, 
s + — ti'P + epx = 


dm 


+ jhvq + jv - \up - \huv + e{qx + jhvy) = 0, 
s + p"^ — wp + j^h?v'^ + epx = 


dm 


pq + ^uq + j+vq^ + e{vqy + qx) = 0, 
s + p^ — v^c^ + wq + e(p2: — t'^gj/) = 


dm 


hpq — ■!(?; + — Y?vc^ + + e/i('i'a; — ^vqy) = 0, 
s + - h?v'^q^ - \h{2w + + e{px - h?v'^qy) = 


dm 


pq - \vq - hvq^ + \u + e{qx - hvqy - (| + hq) Vy) = 0, 

s -p^ + h?v'^q^ + h{huv - hw + - + e{h?v'^{qy - Px) + (1 + 2hq)h?v'^Vx) = 


d3.1()|) 


+ - /i^i) - hvq'^ + X - 1^9(1 - /in) + e {qx - hvqy -h{^ + q) Vy) = 0, 
s+p^ - h?q^v'^ + ^huv"^ + \hq{2hu - l)f ^ - + e (p^ - h'^v'^qy - 2h? + g) wuj/) = 


d3.11j) 


pq - hvq^ - \up + \v + e{qx - hvqy) = 0, 

s+p^ - h^v'^q^ -wp+ \hv'^ + e{px - h^v'^qy) = 



Table 2: Nonlinear pairs for systems p.Sp - ljS.lip . 



4 Conclusion 

The main objective of this paper was to present new classification results based on the method of 
hydrodynamic reductions and their dispersive deformations. The technique was extended to provide 
an algorithm to compute the Lax pair of a dispersive system from that of its dispersionless limit. 
The classification of dispersionless systems was based on the requirement that in order to possess an 
integrable dispersive deformation, a dispersionless system should be generated by a polynomial Lax 
pair. In this work, we have restricted ourselves to second degree polynomial generators and a complete 
list of integrable Davey-Stewartson type systems has been obtained in this context. Consideration of 
higher degree polynomials may also lead to new systems and this task will be undertaken in a future 
work. 
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